Accepted for publication in The Astrophysical Journal 
Preprint typeset using WT^^ style cmulatcapj v. 10/09/06 



THE SPITZER SURVEY OF THE SMALL MAGELLANIC CLOUD: DISCOVERY OF EMBEDDED 

PROTOSTARS IN THE H II REGION NGC 346 

Joshua D. Simon\ Alberto D. Bolatto", Barbara A. Whitney'', Thomas P. Robitaille" , Ronak Y. Shah^, David 
Makovoz'', Snezana Stanimirovic", Rodolfo H. Barba", and Monica Rubio** 

Accepted for publication in The Astrophysical Journal 



o 
o 

(N 



(N 

6 



> 

oo 
a^ 

On 
rn 

l> 
O 

o 



X 



ABSTRACT 

We use Spitzer Space Telescope observations from the Spitzer Survey of the Small Magellanic Cloud 
(S^MC) to study the young stellar content of N66, the largest and brightest H II region in the SMC. 
In addition to large numbers of normal stars, we detect a significant population of bright, red in- 
frared sources that we identify as likely to be young stellar objects (YSOs). We use spectral energy 
distribution (SED) fits to classify objects as ordinary (main sequence or red giant) stars, asymptotic 
giant branch stars, background galaxies, and YSOs. This represents the first large-scale attempt at 
blind source classification based on Spitzer SEDs in another galaxy. We firmly identify at least 61 
YSOs, with another 50 probable YSOs; only one embedded protostar in the SMC was reported in 
the literature prior to the S^MC. We present color selection criteria that can be used to identify a 
relatively clean sample of YSOs with IRAC photometry. Our fitted SEDs indicate that the infrared- 
bright YSOs in N66 have stellar masses ranging from 2 Mq to 17 M©, and that approximately half of 
the objects are Stage II protostars, with the remaining YSOs roughly evenly divided between Stage I 
and Stage HI sources. We find evidence for primordial mass segregation in the H II region, with the 
most massive YSOs being preferentially closer to the center than lower-mass objects. Despite the low 
metallicity and dust content of the SMC, the observable properties of the YSOs appear consistent 
with those in the Milky Way. Although the YSOs are heavily concentrated within the optically bright 
central region of N66, there is ongoing star formation throughout the complex and we place a lower 
limit on the star formation rate of 3.2 x 10^'^ Mq yr^^ over the last ~1 Myr. 

Subject headings: Magellanic Clouds — H II regions — infrared: stars — ISM: individual (N66) — 
stars: formation — stars: pre-main sequence 



1. INTRODUCTION 

Star formation is one of the most important astrophys- 
ical processes, but because it takes place on small physi- 
cal scales and behind heavy optical obscuration, observa- 
tions of the early phases of star formation have generally 
been limited to the Milky Way. Only in recent years 
have these studies begun to be extended to extragalactic 
objects. Now, the high sensitivity, high angular reso- 
lution, and large field of view provided by the Spitzer 
Space Telescope at mid-infrared wavelengths open a new 
window on recently formed stars in nearby galaxies like 
the Magellanic Clouds. Studying these young stellar ob- 
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jects (YSOs) provides an excellent complement to ob- 
servations of star formation in the Milky Way because 
the physical conditions they are subject to are quite 
different from those in the Galaxy. In particular, the 
metallicity of the Magellanic Clouds is well below solar, 
and their dust content is correspondingly l ower as well 
(|Stanimirovic et al.l[2000t iLerov et al]l2007f l. Most star 
formation today is taking place in galaxies smaller than 
the Milky Way (for which the Magellanic Clouds are good 
prototypes), and star formation at high redshift occurred 
in low-metallicity environments, so investigating the ef- 
fects of these differences on how star formation works will 
be an important step towards understanding how many 
of the stars in the universe were formed. 

The Small Magellanic Cloud (SMC), and its more 
massive companion the Large Magellanic Cloud (LMC), 
are the two nearest star- forming dwarf galaxies. 
Loca t ed at a distance of 61.1 kpc dWe sterlundl 
1997; 'Storm et al." '200?; Hilditch. Howarth. fc Harriet 
2005; Keller & Wood 200^, and with a luminos- 
ity of - 6 X 10* Lq 
and a metallicity of ' 



(jde Vaucouleurs et al. 


11991) 


1/5 solar (iDufour 


1975: 



Harlow! 



1977t iPeimbert. Peimbert. fc Rui3 l2000f ). the SMC is 
perhaps the best local analog of primitive galaxies. The 
SMC is actively for ming stars at a rate of ~ 0.05 Mq yr~^ 
()Wilke et al.ll2004f ). and is populated by well studied H II 
regions and young star clusters. The most active star- 
forming region is NGC 346 (alternately known as N66; 
Hcnizc 1956), located toward the northern end of the 
SMC bar. The OB association powering N66 contains 33 



2 



Simon et al. 



spectroscopically confirmed O stars, and a similar num- 
ber of bright blue s tars without spectra that are likely to 
be O stars as well (jMassev. Parker. &: Garmairvlll989D 

In this paper, we employ Spitzer observations from the 
Syitzer Survey of t he Small Magellanic Cloud (S^MC; 
iBolatto et al.ll2007l . hereafter B07) to locate objects in 
N66 with mid-infrared excesses that we identify as can- 
didate YSOs. In the following section, we very briefly 
describe the observations and our photometry. In fj3l we 
classify the detected sources by fitting their spectral en- 
ergy distributions and study their locations in color-color 
space. In S]4] we discuss some of the implications of our 
results and compare to recent optical studies of N66. We 
summarize our findings in !j5] 

2. OBSERVATIONS, DATA REDUCTION, AND 
ANALYSIS 

The S^MC is a project to map the star-forming body of 
the SMC with Spitzer in all seven Infrared Array Cam- 
era (IRAC; Fazio et al. 2004) and Multiband I maging 
Photometer for Spitzer (MIPS: lRieke et al.l[200l bands. 
The images cover an area of ~ 3 deg^ including the entire 
bar and wing of the SMC, with an average exposure time 
of 144 s at each position. The MIPS data were obtained 
in 2004 November and the IRAC data in 2005 May. The 
data were processed with version SI 1.4 of the automated 
Spitzer pipeline. We constructed mosaic images from the 
individual Basic Calibrated Data frames using the Mo- 
saicking and Point Source Extraction (MOPEX) software 
provided by the Spitzer Science Center (SSC). B07 de- 
scribe further details of the observations and the data 
processing. 

We performed photometry on the mosaic images with 
the Astronomical Point Source Extraction (APEX) tasks 
in the MOPEX package (iMakovoz fc Mar lcau 2005). We 
selected a set of 20 — 30 bright stars in each band that 
were as isolated as possible and constructed point re- 
sponse functions (PRFs) directly from the data. We 
then fit these PRFs to every detected source in the im- 
ages to determine fluxes. The images contain exten- 
sive diffuse emission that must be separated from point 
sources, so we used a small median filter (8'.'4) to re- 
move the background, and then detected sources on the 
background-subtracted image. Because the median fil- 
tering removes some flux even from point sources, the 
PSF fitting then took place on the original (not back- 
ground subtracted) image. We found that this technique 
offered the best compromise between detecting bright 
sources on top of diffuse emission and detecting faint 
sources in background-free regions. For the 24 /xm MIPS 
data, in which the extended emission dominates over the 
point sources, we used a somewhat larger median filter 
(27'.'5) and switched to the 'combo' algorithm for sepa- 
rating clusters of bright emission into individual sources 
in place of the 'peak' algorithm in the task detect.^" 
We used observations of bright, isolated point sources 
in the SMC to measure the aperture corrections for the 
photometry. To provide near-infrared fluxes, we used 
VLT/ISAAC J and imaging of the center of N66 

More information about the DETECT tasl< and its associated 
options can be found in tlie documents titled APEX User's Guide 
and Image Segmentation (both by D. Makovoz) that are available 
on the Spitzer Science Center website. 



(|Rubio et al.ll200^ . Rubio & Barba, in p reparation) and 
Two Micron All Sky Survey (2MASS; ISkrutskie et~an 
[20061) measurements in the outer regions. We also added 
HST V- and J-band p hotometry of the central region 
from ^ Gouliermis et al.l (2006), but because the much 
higher angular resolution of the ACS images can cause 
confusion within the Spitzer resolution element, we only 
used these data for bright stars iV < 17.5). 

The assumed photometric uncertainty associated with 
each flux measurement is very important for determining 
the relative weights given to the various data points in 
the SED fitting that we carry out in H3.2[ Although the 
statistical uncertainties on all of the photometric mea- 
surements were quite small (generally a few percent or 
less) , we imposed larger minimum uncertainties on all of 
the data points to account for systematics. For example, 
the HST, 2MASS, VLT, IRAC, and MIPS observations 
were made at various times over a period of 7 years, pro- 
viding plenty of time for the fluxes of the YSOs, which 
may be variable, to chan ge. The absolute photometric 
accuracy of IRAC is 10% dPazio et al.ll200l . so we added 
10% of the flux in quadrature to the measured uncer- 
tainty f or each source. MI PS also has a stated accuracy 
of 10% (jRieke et al.ll2004[ ). so when combined with the 
IRAC calibration accuracy and the 6 month time base- 
line between the observations, we used a 20% minimum 
error for MIPS fluxes. Since the near-IR observations 
took place 6 — 7 years earlier, we imposed larger min- 
imum errors of 35% (for the deep VLT data) and 50% 
(for the shallow 2MASS data) on those data points. The 
optical measurements are more recent (2004 July) and 
more accurate, so we assumed 10% uncertainties on the 
HST photometry. 

3. RESULTS 

3.1. Photometry 

For the purposes of this study, we limited our point- 
source fitting to a 12f75 x 12.'75 region covering the full 
extent of the hot dust emission from N66 seen in our 
24 urn image (from 00''57™5P' to 0l''00™37" and from 
-72°16'48" to -72°04'19" [aU coordinates J2000.0]). In 
this box, we detect 6544 sources at 3.6 /xm, 5836 at 
4.5 /im, 1784 at 5.8 /im, 1718 at 8.0 ^m, 101 at 24 ^m, 
and 15 at 70 fim, for a total of 8011 unique objects. 
A significant number of these objects have Spitzer col- 
ors that are redder than normal stars should be at these 
wavelengths, suggesting that they may be YSOs or back- 
ground galaxies (B07). 

3.2. SED Fitting 

There are two possible approaches to determining the 
nature of individual Spitzer sources. One could simply 
use the observed colors and magnitudes to classify the 
sources into various categories. Stars of nearly all kinds 
(except those with dusty atmospheres) have colors near 
zero for any combination of Spitzer bands because the 
IRAC and MIPS bands are on the Rayleigh-Jeans tail 
of their spectral energy distributions (SEDs).^^ YSOs, 
because of the emission from warm (T ^ 200 K) dust 
around the central protostar, have red colors through- 
out the mid-IR. IRAC colors for theoretical YSO models 

Note that we use the Vega magnitude system throughout this 
paper. 
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are giv en b\^e.g.. lAiren et all (|2004D and lWhitnev et all 
(|2003al lbl. [200l . Alternatively, one can compare the full 
SEDs to a variety of source models and find the best 
match for each object. 

The advantage of the color selection strategy is its sim- 
plicity, but it also has some drawbacks. It fails to use all 
of the available information about each object (since we 
also have MIPS, optical, and near-IR fluxes for many of 
the sources), and it is very difficult to learn about the 
detailed properties of individual sources from only their 
colors. Color selection also does not offer a way to test 
the assumption that SED models designed for Milky Way 
sources offer a good description of YSOs that form in the 
metal- and dust-poor environment of the SMC; this is a 
disadvantage for studies seeking to determine how star 
formation proceeds at low metallicity, but could also be 
an advantage in that it enables the selection of YSOs 
without regard to the properties of their environment. 
Because the SED classification should be more accurate, 
we begin with that technique and then compare the re- 
sults to color selection in i i3.3l 

The SED fitting tool used for this study employs a 
linear regression method to find all the SEDs from a 
large grid of models that fit th e data within a spec- 
ified (|Robitaille et al.l I2007D . The grid of mod- 
els consists of 7853 stella r atmospheres (iKurucd Il992l : 
iBrott fc HauschildtJ [2005h encompassing all available 
metallicities and effective temperatures, a limited num- 
ber of Infrared Space _ Observatory (ISO) spectra of galax- 
'Silva et al) "igOS") and AGB stars 
1199 9; Olivier et al. 2001; Honv 2002; 



10" 



ies ("Dale et aTl 



([S ylvester et al., ^^^^ 
Molster et al.l l2002t illonv et all l2002l : iFuiii fc Nakada 



20031) and the 20,000 YSO models from lRobitaille et al 



(l2006t) computed using the rad iation transfer codes from 
iWhitnev et all (|2003al lbl. [200l . Each YSO model out- 
puts SEDs for 10 viewing angles, so the YSO grid ef- 
fectively contains 2 x 10'"' SEDs. The foreground ex- 
tinction, Ay, is fit simultaneously using an extinction 
law derived from GLIMPSE observations (Indebetouw 
et al. 2006). At ultraviolet and visible wavelengths, 
this extinction law is not appropriate for the SMC, but 
in the near-IR and mid-IR the differences in extinction 
between the SMC and Gal axy are small (jGordon et al.l 
[20031: iCartledge et al.l[2?)05h . The fitter is run first using 
only the stellar atmosphere grid. It is then run three 
more times using the YSO grid, ISO galaxy spectra, and 
ISO AGB spectra on the sources that are not well-fit 
by stellar atmospheres. Based on all the successful fit 
results, defined by their values (see below), we can 
classify sources and calculate best estimates and uncer- 
tainties for each model parameter. 

We fit an SED for every object for which we had at 
least four flux measurements. To better constrain the 
fits at long wavelengths where most of the sources were 
not detected, we added 24 /im upper limits of 1 mJy (5 
times the limiting sensitivity of the 24 /im data) for each 
source that was not detected at 24 /tm. Because there 
are many more model parameters than data points that 
go into each observed SED, the traditional reduced 
statistic cannot be used to evaluate the success of the fits. 
As a proxy, we calculate the per data point (where 
each SED has between 4 and 11 data points). Given this 
definition, it is difficult to establish a priori the division 
between good and poor fits. We inspected the fits by 
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Fig. 1.— Sample SED fits for a YSO in N66. The black points 
indicate the measured fluxes and uncertainties, and the black tri- 
angles are for upper limits. The thin black line represents the 
best fitting SED, and the gray lines represent all other acceptable 
(x^/data point < 2.2) YSO fits. The dashed line represents the 
stellar photosphere model (including the effect of foreground ex- 
tinction). T his o bject is a Stage II YSO with a stellar mass of 
4.9 M0 (see for the definition of YSO stages). 

eye and determined that fits with a x^ psr data point 
of less than 2.2 could reasonably be interpreted as good 
fits, while higher x^ values indicated fits that were not a 
good match to the templates. We will therefore use this 
cutoff point to separate good and poor fits throughout 
the paper. Note that this definition depends on the fiux 
uncertainties imposed in ij2l 

3.2.1. Results From Successful SED Fits 

Out of the 1645 sources with four or more ffuxes, 1322 
(80.4%) are succesfully fit with stellar SEDs (see above 
for what we mean by a "successful fit"). Of the remain- 
ing 323 sources, 61 can only be fit successfully by YSO 
models, 27 are background galaxies, 6 are AGB stars, 81 
can be fit by multiple classes of models, and 148 cannot 
be fit by any of the available templates. We display one 
example YSO SED in Figure ffl 

For the sources that can be well-described by multiple 
types of templates, we used the x^ value of the best fit for 
each template type to attempt a more accurate classifica- 
tion. We set a cutoff of 1.5 times the best x^ value, and if, 
for example, the best YSO fit had x^ = 1-0 and the best 
galaxy fit had x^ = 1-6 then we classified the object as a 
probable YSO. This choice of a x^ cutoff indeed results in 
the majority of the uncertain objects being classified as 
YSOs and leaves very few objects with unknown types, 
which appears to be a reasonable outcome. While the 
exact value chosen for the cutoff is arbitrary, it is clear 
from the spatial distribution of the various object types 
that many of the sources classified as possible or defi- 
nite galaxies must actually be YSOs. The results of this 
paper do not depend strongly on the exact value of the 
cutoff. Using this method, we were able to obtain prob- 
able classifications for 68 of the 81 objects of uncertain 
type, including 50 YSOs, 13 galaxies, and 5 AGB stars. 
Thus, the final object counts are 111 YSOs, 40 galax- 
ies, and 11 AGB stars. The remaining 13 objects had 
X^ values for multiple object types that were too close 
to distinguish reliably. We list selected parameters from 
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the fits for all 111 YSOs in Table [H and we summarize 
the source classification results in Table [TJ 

3.2.2. Objects Poorly Fit by the SED Models 

We inspected each of the 148 poorly fit sources and the 
various attempted fits to their SEDs to determine why 
the fitting failed for them. We found that slightly more 
than 1/3 of the objects had SEDs at short wavelengths 
that appeared stellar in origin, but the long wavelength 
data (often at 24 /im) were significantly in excess of the 
extrapolated photospheric emission. The fits to these 
sources appear to have failed for a number of reasons 
including slight mismatches between the optical and IR 
photometry (either from variability or measurement er- 
ror), confusion, misidentification of the long wavelength 
counterparts, and contamination of the photometry by 
underlying dust emission (revealed by PAH features). 
However, some of these sources may also have true in- 
frared excesses indicating the presence of circumstellar 
material. Another 1/3 of the bad fits were gener- 
ally faint sources that are spatially coincident with dust 
filaments in the Spitzer images, again causing signifi- 
cant contamination of their SEDs by PAH features in 
the IRAC bands. Some of these objects are probably 
point sources (either stars or YSOs), while others may 
simply be unresolved dust knots. The final 1/3 of the 
bad fits was a mixed population whose SEDs could not 
be straightforwardly interpreted. Photometric errors in 
one or more bands and blending were likely responsible 
for the failure to fit these objects, but they represented 
only 2.6% of the total sample. 

3.3. IRAC Colors of YSOs 

Using the results of the SED fitting, we can now in- 
vestigate the colors of the objects classified as YSOs and 
as stars. We display four color-color plots in Figure [2] 
to illustrate the possibilities for color selection. Stars, as 
expected, lie in the cloud of points centered near (0, 0), 
while sources that are red in one or both colors may 
be YSOs. Note that the frequently used [3.6] -[4.5] vs. 
[5.8] — [8.0] color-color plot (Fig. ^) does not cleanly 
separate YSOs from other types of sources, particularly 
stars with modest IR excesses and sources with PAH 
contamination. Better separation can be achieved using 
different combinations of colors, as shown in Figure ^p. 
These plots take advantage of longer color baselines and 
the abrupt change in YSO spectra between the 4.5 fim 
and 5.8 /im bands to distinguish YSOs from stars, galax- 
ies, and PAHs. These results demonstrate that while 
color selection can be a useful technique for identifying 
YSOs, it does not appear possible to obtain a YSO sam- 
ple that is both complete and clean with a simple set of 
color diagnostics. SED fitting is a more comprehensive 
way to determine the nature of sources if measurements 
in enough bands are available. Nevertheless, because in 
many cases it is desirable to classify objects with easily 
applicable techniques, we use our data to provide guid- 
ance for color selection of YSOs. 

Assuming that the purple points in Figure [2] are in- 
deed YSOs, one can select YSOs with the following set 
of criteria: 

[3.6] - [4.5] > 0.6 X ([4.5] - [8.0]) - 1.0 



[4.5] - [8.0] < 2.8 

[3.6] - [4.5] < 0.6 X ([4.5] - [8.0]) + 0.3 

[3.6] - [4.5] > -([4.5] - [8.0]) + 0.85. (1) 

These criteria were defined to maximize completeness; 
one could alternatively choose to minimize contamina- 
tion at the cost of increased incompleteness, but given 
the distribution of sources in Figure the differences 
would be small. This selection is only 7% incomplete for 
YSOs in our sample that have measured fluxes at 3.6 fim, 
4.5 /im, and 8.0 /im, and has a contamination of less than 
27% (the majority of the objects in this color box with 
uncertain classifications may still be YSOs). Approxi- 
mately equivalent results can be obtained by substitut- 
ing the [3.6] -[8.0] color for the [4.5] -[8.0] color, but it 
is slightly more difficult to avoid picking up stars on the 
blue end of the [3.6] — [8.0] color axis. 

3.3.1. Comparison to B 07 Photometric Selection 

B07 identified a set of very bright YSOs across the 
entire SMC based on their 8.0 /xm magnitudes and 
[5.8] — [8.0] colors. If we select the same region of the 
color-magnitude diagram in N66, we find 15 sources, only 
3 of which are conclusively classified as YSOs. Almost all 
of the remaining objects have SEDs that are not well-fit 
by any of the models we apply, primarily because of PAH 
contamination. The PAH contamination across the rest 
of the galaxy should be much lower than in N66, so these 
numbers clearly represent a lower limit to the fraction of 
actual YSOs in the B07 sample. 

3.4. YSO Masses and Stages 

YSOs are traditionally divided into classes based on 
their observ ed spectral indices as originally defined by 
iLadal ()1987n . Because spectral indices can vary with 
inclination angl e as w ell as with evolutionary state, 
iRobitaille eFall (|2006D described their YSO models in 
terms of "stages", which are analogous to the usual 
classes, but are based on the physical quantities that 
define the evolutionary stage of the models. When 
comparing these models to data (rather than consider- 
ing only observations), it therefore makes sense to use 
the stage system. The definitions of each stage rely 
on the ratio of the disk mass and envelope accretion 
rate to the central stellar mass. Stage I sources have 
Af/M, > 10-6, Stage II sources have M/M^ < lO'^ 
and Mdisk/M^ > 10^^, and Stage HI sources have 
M/M, < 10-6 and Mdisk/M^ < 10'^, where M is the 
envelope accretion rate. Note that we assumed a dust- 
to-ga s ratio of 2 x 10"^ (1/5 of the Milky Way value) for 
N66 (iLerov et al. 1 12001 .1^ 

Using this classification scheme, we have grouped the 
111 definite and probable YSOs in N66 into stages. The 
SEDs of most objects can be fit by more than 1 YSO 
model, so we first selected the models that produce a 
that is within 1 of the best for each object. We then 
computed a stage for each of the acceptable models and 
calculated a weighted average stage, using the values 
as weights. The averages were rounded to the nearest 
integer to produce a classification. We find that of the 

The dust-to-gas ratio does not affect the YSO fits themselves, 
but does change the mass ratios between the dislt and the central 
star, which is used for the stage classification. 
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Fig. 2. — (a) IRAC [3.6] — [4.5] vs. [5.8] — [8.0] color-color diagram of N66. The black dots represent all detected sources (primarily 
stars), the filled red circles represent objects that can only be fit by YSO models, the smaller purple circles represent objects that are best 
fit by YSO models, although galaxy and/or AGB fits with significantly higher values are also acceptable, the filled blue triangles are 
background galaxies, and the filled green diamonds are AGB stars. The open squares represent sources for which satisfactory fits were not 
obtained: the cyan symbols have SEDs that suggest they are likely to be stars, the orange symbols are sources that are contaminated by 
PAH emission features, and the gray symbols are the remaining unclassified poor fits. Note that the separation of YSOs from the other 
classes of sources is not very clean in this diagram. (6) IRAC [3.6] — [4.5] vs. [4.5] — [8.0] color-color diagram of N66. Symbols are the same 
as in (a). This plot offers the best separation of YSOs from the various contaminants, although there are still a few extragalactic sources 
that have similar colors to the YSOs. 



Ill YSOs, 33 are Stage I, 50 are Stage II, and 28 are 
Stage III. 

We also calculated YSO masses with the same weighted 
averaging scheme, and defined the uncertainty on the 
mass to be the weighted standard deviation of the masses 
of the acceptable models. We found YSOs with a range 
of masses from 2.4 — 16.6 Mq, including 19 objects that 
appear to be proto-OB stars (M > 8 M0). We list the 
names, positions, luminosities, luminosity uncertainties 
(AL, the weighted standard deviation of the luminosi- 
ties of the acceptable models), masses, mass uncertainties 
{AM, as defined in the first sentence of this paragraph), 
stages, and stage uncertainties (AStage, the weighted 
standard deviation of the stages of the acceptable mod- 
els) of all of the YSOs in Table □ 

3.5. Spatial Distribution of YSO Candidates 

In Figure[3]we plot the spatial distribution of the YSOs 
in N66. The protostars are obviously highly concentrated 
towards the peaks of the 8.0 yum emission. However, there 
is also star formation taking place outside of the dense 
dust cloud that marks the optical H II region, particu - 
larly to the south and southeast. llVIassev et al.l ()1989D 
hypothesized that star formation in N66 has proceeded 
from the southwest to the center of the present-day H II 
region, but we do find a few probable YSOs southwest of 
N66, indicating that at least modest star formation has 
taken place there within the past few million years. 

The distribution of YSOs throughout the H II region 
as a function of mass and stage is not uniform. We find 
that the most-embedded objects (Stage I) are slightly 
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Fig. 3. — Spatial distribution of YSOs in N66, overplotted on an 
image of the 8.0 ^tm emission. The red circles represent objects with 
SEDs that can only be fit by YSO models, and the yellow circles 
represent objects that are probably YSOs but can be fit by other 
source types as well (albeit with significantly higher values). 
The image is displayed on a logarithmic scale from 4 MJy sr~^ to 
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Fig. 4. — Mass distribution of YSOs in N66, overplotted on an 
image of the 8.0 /im emission. The most massive sources (M > 
8 Mq) are plotted as cyan points, and then decreasing masses are 
indicated by yellow (4.5 M0< M < 8 Mq) and red (M < 4.5 Mq) 
points. The high-mass stars are more concentrated in the center 
of the H II region than the less massive objects. Note that the 
objects shown in this figure include the full YSO sample, not just 
the 61 definite YSOs that can only be fit by YSO SEDs. 

more concentrated towards the center of the H II region 
than the more advanced (and presumably older) YSOs. 
We also see evidence for mass segregation, with the most 
massive objects exhibiting a strong preference for loca- 
tions close to the center (see Figure H]) . All but 2 of the 
YSOs with M > 8 M© he on top of bright dust fila- 
ments in the main H II region, and many of them arc 
coin cident with molecular peaks and op tical star clus- 
ters (|Rubio et al.ll2000t ISabbi et al.ll2007D . If these mas- 
sive YSOs are indeed single objects rather than multiple 
unresolved sources, then the mass segregation must be 
primordial in origin, as the YSOs have not had time to 
move very far from their birthplaces. 

4. DISCUSSION 

The census of previously known embedded proto- 
stars in the SMC c o nsists of a single object di scov- 
ered bv iGatley et all (|1982D . iBeaulieu et"an (|2001[) and 
Ide Wit et al.l ( 2003f ) used variability data and Ha imag- 
ing to uncover a small samp le of Herbig Ae/Be stars 
in the SMC, and iNota et alj |2006D discovered several 
hundred low-mass pre-main-sequence stars in N66 via 
isochrone fits to HST photometry. These objects, how- 
ever, are generally much more evolved than the YSOs 
identified by this study, many of which are still embedded 
in their natal dust clouds and are therefore faint or in- 
visible at optical wavelengths. The distribution of YSOs 
is concentrated towards clum ps of molecular ga s traced 
by the C0(2-l) emission line (iRubio et al.ll200C)D and th e 
peaks of the dust em ission at 7 /im (jContursi et aL|[200l . 
IRubio et al.l ([2000) showed that dense H2 knots are as- 
sociated with these molecular clumps and suggested that 
massive star formation could be taking place there. The 



large number of embedded YSOs found in this study con- 
firms that prediction. Moreover, IR spectroscopy of the 3 
brightest embedded sources detected with ground-based 
near-IR imaging confirm that these sources are YSOs 
(Rubio and Barba, in preparation). 

The sample of Stage I, II, and III YSOs that we have 
identified in N66 presents the first opportunity for study- 
ing in detail a large sample of e mbedde d YSOs in an- 
other galaxy (see also lJones et al. i 2005; C hu et an [2065). 
and more importantly, one whose ISM properties differ 
substantially from those of the Milky Way. The fact 
that we are able to identify over 100 YSOs in N66, as 
well as obtain successful SED fits for almost all of the 
other sources, suggests (perhaps surprisingly) that proto- 
sta rs in the SMC r esemble the YSO models constructed 
by IWhitnev et al.l for Milky Way objects, even though 
the metallicity and dust-to-gas ratio are a factor of ~ 5 
lower in the SMC. The only potential difference between 
Milky Way and SMC YSOs that is evident in our re- 
sults is that some of the sources are best fit by relatively 
cool photospheres, but the YSO models with such pho- 
tospheres do not have high enough luminosities to match 
their observed brightnesses at the distance of the SMC. 
If these objects are actually YSOs, this problem could 
be an indication that accretion is continuing even af- 
ter these stars have reached the main sequence, which 
leads to expanded photospheres and lower temperatures 
compared to normal pre-main seq uence tracks and zero - 
age main-sequence photospheres (|McKee fc TanI [20031 ) . 
An alternative possibility is that these sources actually 
consist of multiple cool protostars, which would explain 
their unusually high luminosities. The luminosities of 
these objects are 2 — 3 orders of magnitude higher than 
the expected luminosities of individual YSOs of the same 
temperature, however, which makes the multiplicity ex- 
planation appear unlikely. Finally, evolved stars in the 
post-AGB phase also have expanded cool photospheres, 
so it is possible that some of the "cool luminous" sources 
could be post-AGB stars. Additional modeling and ob- 
servations of some representative objects may be neces- 
sary to resolve this issue. 

4.1. The Star Formation Rate in N66 

The 111 YSOs in our sample have a combined stel- 
lar mass of 692 M©. The observed mass function (see 
Figure El) turns over at ~ 4 Mq, indicating that incom- 
pleteness becomes serious at this point. Artificial star 
tests in our photometry show that we are 90% complete 
even in the confused central region of N66 down to flux 
levels of 300 /iJy, 200^iiy, 200/iJy, and 300 ^iJy from 3.6- 
8.0 //m. Applying these limits to the full library of YSO 
models confirms that our incompleteness is severe below 
4 Mq. If we assu me that star formation in N66 follows 
a ISalpeteij (|1955f ) initial mass function (IMF) down to 
0.1 Mq, we calculate that the total mass in protostars 
for the entire H II region is ~ 3160 Mq. In reality, this is 
a lower limit to the mass because even at 4 Mq the data 
are somewhat incomplete. If these YSOs have all formed 
within the last ~ 1 Myr, then the average star formation 
rate over that time is 3.2 x 10"'^ Mq yr~^. Thus, we find 
that N66 comprises at least ~ 6% of the total current 
star formation in the SMC. 

4.2. Comparison with Optical Observations 
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iNota et alJ (l2006l) . iGouliermis et al.l (|2006D . and 
iSabbi et all Mm used HST imaging to stu dy the young 
stellar population of N66 in the optical. ISabbi etldl 
(|2007l ) identified 16 subclusters of pre-main sequence 
stars (their Fi gure 8; several of these clusters were also 
pointed out bv IGouliermis et al.l ) in the H II region, and 
we find that all but two of t hese also hav e YSOs associ- 
ated with them (see below). ISabbi et"al1 estimated ages 
of ^ 3 Myr for the first 15 of the subclusters, but our 
detection of YSOs demonstrates that star formation has 
continued until the present day in these areas. 

Subcluster 1 (Sc 1) is nearly coincident with an ex- 
tremely bright Spitzer source (the 4th-most luminous ob- 
ject in the field at 8.0 nm), SSTS3MC 14.7725-72.1766, 
which we are able to fit with both YSO and AGB SEDs 
with similar values. Given the location of this ob- 
ject at the center of the NGC 346 cluster, surrounded by 
numerous very young massive stars, it is most likely a 
YSO (nevertheless, since it does not formally meet our 
selection criteria it is not included in our analysis). If 
so, the fitted SEDs suggest that this is a Stage I ob- 
ject with a luminosity of 3.3 x 10* Lq and a mass of 
14.7 Mq. However, it is important to remember that 
because of the very high source density here we may ac- 
tually be seeing multiple unresolved YSOs. In that case, 
we would be likely to overestimate the mass of the most 
massive YSO and underestimate the total mass of YSOs 
contained in this source. This type of source confusion 
should not have a strong effect on the SED fitting results, 
because the observed SED will be dominated by that of 
the most massive embedded YSO. A mid-IR spe ctrum 
of this object is displayed in IContursi et al.l (|2000l . peak 
C). Another bright Spitzer source, SSTS3MC 14.7748- 
72.1749, is located within 3" of Sc 1 and has a steeply- 
rising mid- and far-IR SED that we are unable to fit suc- 
cessfully. If we remove the 24 fim upper limit and either 
the 70 fim or J- and K-hand detections (perhaps justi- 
fied because of confusion in this very densely populated 
region), then this source has the SED of a very massive 



(M > 10 Mq), early-stage YSO. We do not detect any 
non-stellar sources in Sc 3 (immediately south of Sc 1), 
but confusion as a result of the very bright source just 
north of the cluster may play a role in this non-detection. 

Sc 2 is located very close to 2 Spitzer sources. 
SSTS3MC 14.7574 has a stellar SED with a strong 24 ^im. 
detection (indicating either confusion or circumstellar 
dust) and SSTS3MC 14.7580-72.1763 is formally a poor 
fit as a result of a low if-band fiux and slight PAH con- 
tamination, but appears to have an SED consistent with 
being an early stage YSO. 

Scs 4-6 each have bright infrared counterparts. Sc 4 
is coincident with SSTS3MC 14.7605-72.1687, a mas- 
sive embedded YSO. Sc 5 is SSTS3MC 14.7514-72.1681, 
another object that can be fitted by both YSO and 
AGB SEDs. Again, given its position in a very 
young cluster, it is most likely a massive stage I YSO 
rather than an evolved star. And Sc 6 contains 2 
blended Spitzer sources, SSTS3MC 14.7371-72.1651 and 
SSTS3MC 14.7380-72.1651, both very high-mass stage I 
YSOs. The caveat mentioned above about multiple un- 
resolved sources applies here as well. 

Located on the southern edge of Sc 7, very close to 
the prominent dust lane that arcs nearly halfway around 
N66, is SSTS3MC 14.7733-72.1835, yet another very 
massive (18.6 Mq) stage I protostar that can also be 
fit by AGB models. Slightly further to the southwest 
are 2 additional, lower-mass YSOs, SSTS3MC 14.7673- 
72.1834 and SSTS3MC 14.7698-72.1846. 

Sc 8 has a faint Spitzer counterpart, 
SSTS3MC 14.7816-72.1802, which shows a stellar 
SED in the optical and near-IR and excesses in the 
IRAC bands, but is not fit well by any of our YSO 
models. It is possible that the optical and IR emission 
is coming from different sources. 

Scs 9-11 are located around a quintet of YSOs. The 
2 brightest of these, SSTS3MC 14.8112-72.1843 and 
SSTS3MC 14.8130-72.1840, lie within Sc 10, but it is 
not clear whether these are truly blended YSOs as 
opposed to a bright, extended clump of dust. The 
other 3 sources classified as YSOs in this region are 
SSTS3MC 14.8041-72.1867, SSTS3MC 14.8068-72.1856, 
and SSTS3MC 14.8174-72.1866. 

Sc 12 contains two bright Spitzer sources 
(SSTS3MC 14.8318-72.1890 and SSTS3MC 14.8353- 
72.1892) that are not fit well by any of the available 
models. The rising SEDs of these objects towards 
long wavelengths strongly suggest that there are YSOs 
present here, but blending may be a problem. This 
cluster is coincident with the source labeled peak I by 
IContursi et all (|2000D . 

Scs 13-15 lie in the northward extension of the H II 
region (N66A). Clusters 13 and 14 each overlap w ith a 

peak 
peak 



Contursi et al 



Contursi et al 



YSO (SSTS3MC 14.8009-72.1663: 
F and SSTS3MC 14.8080-72.1578: 
G), but as with Sc 10 the source in Sc 14 does not appear 
pointlike in the IRAC bands, and so may not be a sin- 
gle object. Sc 15 is associated with SSTS3MC 14.8205- 
72.1544, which shows very strong PAH emission but is 
not well-fit by YSO models. 

Finally, Sc 16 contains 2 Spitzer sources, but both of 
these are well-fit by normal stellar models, consistent 
with t he older age of this cluster derived by iSabbi et all 
(|200l . 
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Out of the 44 YS Os in our sample t hat he within the 
ACS Ha images of iNota et all (|2006[ ). ~ 10 are spa- 
tially coincident with gas or dust pillars strongly resem- 
bling those made famous by HST im aging of the "Pillar s 
of Creation" in the Eagle Nebula (Hes ter et al.l 11996). 
Higher spatial resolution near-IR and mid-IR imaging of 
these objects may reveal exactly how the YSOs are re- 
lated to these features. 

5. SUMMARY AND CONCLUSIONS 

We have obtained mid-infrared imaging of the H II re- 
gion N66 (NGC 346) in the SMC with the IRAC and 
MIPS instruments on the Spitzer Space Telescope. We 
detected 8011 unique sources, with photometric cover- 
age extending from V^-band to 24 /im (70 /im in a few 
cases) . Most of these sources have colors and SEDs con- 
sistent with being normal stars, but we also detect a 
significant population of objects that are very red in the 
mid-IR. SED fitting of the 1645 sources with photometric 
measurements in at least four bands yielded 111 objects 
with SEDs that are best fit by YSO models rather than 
stars or background galaxies. These data represent the 
first significant sample of embedded YSOs identified in 
an external galaxy. 

We show that these YSOs can be mostly, but not com- 
pletely, separated from stars on the basis of their IRAC 
colors. However, SED fitting is necessary to significantly 
constrain the properties of ind ividual o bjects. We find 
that the YSO models of Whitne£eLalJ (2003a,b, 2004), 
which were designed to represent Milky Way YSOs with 
solar metallicities, fit most of the YSOs in N66 well, de- 
spite its much lower metallicity and dust-to-gas ratio. 
These results suggest that if low metallicity causes signif- 
icant changes in the star formation process, the threshold 
for those effects must li e below the me tallicity of the SMC 
(12 + log[0/H] « 8.0: iDufoudflQTSh . The one possible 
difference between SMC and Milky Way YSOs is that the 
SMC appears to contain a population of very luminous 
but cool objects that have not been seen in the Milky 
Way. These objects may have expanded photospheres 
and hence lower temperatures than normal because they 
are still accreting material from their protostellar disks 
when they reach the main sequence. Alternatively, these 
sources could be multiple cool YSOs that are unresolved 
by Spitzer. 

We calculate masses and stages (analogous to the usual 
YSO classes) for each of the N66 YSOs, finding a range 
of masses from 2.4 to 16.6 Mq (including 19 objects with 
masses above 8 Mq). Almost half (45%) of the YSOs are 
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Young Stellar Objects in N66 
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00 59 18.82 


-72 09 24.3 
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5.5 




II 




4.88 




43 
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00 59 19.25 
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III 
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0.00 




44 


SSTS3MC14.8373-72.1449 


00 59 20.95 


-72 08 41.5 
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345 
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0.4 


II 


0.3 


0.15 
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SSTS3MC14.8465-72.1468 


00 59 23.16 


-72 08 48.4 


4430 
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III 
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2.67 




46 


SSTS3MC14.8487-72.1464 


00 59 23.69 


-72 08 47.1 


204 


60 


3.9 


0.3 


II 




2.78 




47 


SSTS3MC14.8551-72.1345 


00 59 25.22 


-72 08 04.1 


1340 




6.6 




III 




0.40 




48 


SSTS3MC14.8701-72.2089 


00 59 28.83 


-72 12 31.9 


131 




5.1 




I 




6.78 




49 


SSTS3MC14.8830-72.2265 


00 59 31.91 


-72 13 35.3 


206 




5.9 




I 




5.28 




50 


SSTS3MC14.8900-72.2789 


00 59 33.61 


-72 16 43.9 


678 


215 


5.5 


0.5 


III 




2.11 




51 


SSTS3MC14.9039-72.1509 


00 59 36.94 


-72 09 03.2 


77 


11 


4.6 


0.3 


I 


0.4 


3.33 





TABLE 1 — Continued 



Number 


Source name 


a (J2000.0) 


S (J2000.0) 


L(Lo) 


AL- (Lq) 


M. (Mq) 


AM*^ (M©) 


Stage 


AStage^ 


^2b 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


52 


SSTS3MC14. 9041 — 72.2398 


00 59 36.98 


— 72 14 23.2 


1000 


495 


5.8 


1.5 


II 


0.2 


0.22 




53 


SSTS3MC14. 9360— 72.2122 


00 59 44.64 


— 72 12 44.0 


2150 




7.5 




III 




7.94 




54 


SSTS3MC14.9603-72.2764 


00 59 50.48 


-72 16 35.2 


644 


249 


5.4 


0.6 


III 


0.5 


0.51 


— 


55 


SSTS3MC14.9671-72.2613 


00 59 52.12 


-72 15 40.8 


529 


376 


5.0 


0.9 


II 


0.5 


0.21 




56 


SSTS3MC 14.9858- 72. 2454 


00 59 56.59 


-72 14 43.5 


683 


177 


5.9 


0.8 


II 


0.4 


11.79 




57 


SSTS3MC15.0636-72.2746 


01 00 15.27 


-72 16 28.7 


1000 




6.2 




III 




5.99 




58 


SSTS3MC15.0655-72.0788 


01 00 15.72 


-72 04 43.6 


360 


87 


4.6 


0.4 


II 




0.26 




59 


SSTS3MC15.0859-72.2232 


01 00 20.62 


-72 13 23.6 


488 


119 


5.0 


0.3 


II 


0.5 


7.91 




60 


SSTS3MC15.0930-72.1662 


01 00 22.32 


-72 09 58.2 


2910 


944 


8.0 


0.7 


I 


0.5 


0.40 




61 


SSTS3MC15.1415-72.2574 


01 00 33.97 


-72 15 26.5 


1750 


99 


8.8 


0.4 


I 




5.55 
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92 


SSTS3MC14.7893-72.1401 


00 59 09.43 


-72 08 24.2 


80 


40 


3.6 


0.5 


II 


0.5 


0.10 


7.14 


93 


SSTS3MC14.7900-72.0911 


00 59 09.60 


-72 05 27.9 


821 


226 


5.8 


0.4 


III 


0.7 


0.48 


3.25 


94 


SSTS3MC14.7908-72.2089 


00 59 09.80 


-72 12 31.9 


57 


23 


2.7 


0.2 


II 




0.38 


8.03 


95 


SSTS3MC14.7909-72.1970 


00 59 09.81 


-72 11 49.1 


2270 


1610 


7.2 


1.7 


III 


0.5 


0.49 


9.59 


96 


SSTS3MC14.7911-72.2052 


00 59 09.87 


-72 12 18.8 


162 


74 


3.7 


0.3 


II 


0.3 


0.05 


6.49 


97 


SSTS3MC14.7956-72.1865 


00 59 10.93 


-72 11 11.5 


114 


373 


3.2 


1.5 


II 




0.01 


3.12 


98 


SSTS3MC14.7988-72.1372 


00 59 11.72 


-72 08 13.7 


114 




5.0 




I 




6.32 


12.05 


99 


SSTS3MC14.8006-72.1939 


00 59 12.13 


-72 11 38.1 


171 


54 


4.1 


0.7 


I 


0.5 


1.69 


8.27 


100 


SSTS3MC14.8011-72.1935 


00 59 12.25 


-72 11 36.7 


305 


228 


4.2 


0.8 


II 




0.07 


7.34 


101 


SSTS3MC14.8426-72.2461 


00 59 22.23 


-72 14 45.8 


275 


170 


4.7 


0.6 


II 


0.2 


0.00 


0.12 


102 


SSTS3MC14.8466- 72.1513 


00 59 23.19 


-72 09 04.5 


356 


109 


4.5 


0.6 


II 




2.22 


6.76 


103 


SSTS3MC14.8532-72.1669 


00 59 24.77 


-72 10 00.8 


1160 


90 


7.0 


0.7 


II 


1.0 


0.70 


7.04 



TABLE 1 — Continued 



Number 


Source name 


a (J2000.0) 


S (J2000.0) 


L(Lo) 


AL- (Lq) 


M. (Mq) 


AM,^ (Mq) 


Stage 


AStage^ 


^2b 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


104 


SSTS3MC14. 8715— 72.1774 


00 59 29.15 


— 72 10 38.6 


490 


58 


7.3 


0.3 


I 




2.46 


5.53 


105 


SSTS3MC14.8729-72.1393 


00 59 29.49 


-72 08 21.6 


123 


45 


4.4 


0.8 


I 




2.10 


4.92 


106 


SSTS3MC14.8781 -72. 1694 


00 59 30.74 


-72 10 09.9 


441 


385 


4.4 


1.5 


III 


0.5 


0.97 


9.78 


107 


SSTS3MC14.8817-72.2385 


00 59 31.60 


-72 14 18.6 


540 


963 


4.6 


1.6 


II 


0.7 


0.16 


6.91 


108 


SSTS3MC14.8848-72.0817 


00 59 32.35 


-72 04 54.2 


98 


47 


4.2 


1.1 


I 




1.85 


6.99 


109 


SSTS3MC14.9738-72. 1101 


00 59 53.70 


-72 06 36.3 


201 


215 


4.3 


0.7 


II 


0.6 


0.05 


4.07 


110 


SSTS3MC15.0896-72.1743 


01 00 21.49 


-72 10 27.5 


80 


30 


3.9 


0.9 


I 


0.5 


0.14 


3.10 


111 


SSTS3MC15.1199-72.1377 


01 00 28.78 


-72 08 15.6 


847 


166 


5.9 


0.3 


III 




0.06 


6.27 



Note. — The coordinates in this table are not necessarily identical to the ones presented by B07 and made available on the S^MC website because the N66 photometry was carried out 
independently. However, the differences should be very small. 

^ The uncertainties on the luminosities, masses, and stages are calculated as the weighted standard deviation of the luminosities, masses, and stages of all of the acceptable YSO models. In cases 
where there is only one acceptable model for a given source, or there are multiple models but they all produce the same luminosity/mass/stage, we cannot calculate an uncertainty in this way, so 

we leave the corresponding space in the table blank. ^"^ x"^ value for the best fitting SED. Note thai these are raw values, not per data point, so a good fit {x^ per data point < 2.2) can have 
a x^ high as 8.8 (if there are 4 flux data points for that source) or more.^ x^ value for the best fltting AGB star or background galaxy SED. 
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TABLE 2 
SED Classification Results 



Object Class Number 

All sources 1645 

Objects with unique SED classifications 

Normal stars 1322 

YSOs 61 

Galaxies 27 

AGB stars 6 

Objects without unique SED classifications 

Probable YSOs 50 

Probable galaxies 13 

Probable AGB stars 5 

Multiple good fits; unable to classify 13 

Poor fits 148 



